Quench Dynamics of Edge States in 2-D Topological Insulator Ribbons 
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We study the dynamics of edge states of the two dimensional BHZ Hamiltonian in a ribbon ge- 
ometry following a sudden quench to the quantum critical point separating the topological insulator 
phase from the trivial insulator phase. The effective edge state Hamiltonian is a collection of decou- 
pled qubit-like two-level systems which get coupled to bulk states following the quench. We notice 
a pronounced collapse and revival of the Loschmidt echo for low-energy edge states illustrating the 
oscillation of the state between the two edges. We also observe a similar collapse and revival in the 
spin Hall current carried by these edge states, leading to a persistence of its time-averaged value. 



PACS numbers: 



:.70.Tg, 03.65.Pm, 74.40.Kb 



Topological insulators (TIs) are novel materials with 
an insulating bulk and conducting edges which are of 
extensive contemporary interest [IHl]. The low-energy 
electrons in two dimensional (2-D) Hg-Te/Cd-Te quan- 
tum well TIs, which display conducting helical edge state 
solutions that exist within the bulk bandgap in the TI 
phase, are described by the 2-D BHZ Hamiltonian [51 H] . 
The bulk states undergo a quantum phase transition 
(QPT) [SHi] with the low-energy modes satisfying a 2- 
D Dirac Hamiltonian (DH) with a linear dispersion at 
the quantum critical point (QCP) (which is a 2-D Dirac 
point). Additionally, the chiral edge states with linear 
dispersion in the TI phase are described by an effective 

1-D DH [g. 

At the same time, there is a recent upsurge in stud- 
ies of quenching dynamics of quantum many body sys- 
tems across QCPs [BHTO]. essentially because of the pos- 
sibility of experimental realization of the same in optical 
lattices [11]. The scaling of the defect density gener- 
ated in the final state following a slow [TH [T3] or a sud- 
den quench |14l I15j , or generation of quantum correla- 
tions which are otherwise absent in the defect free final 
state [TB] or the possibility thermalization with an effec- 
tive temperature I17j are some of the topics which are 
being explored thoroughly. 

In this communication, we focus on the dynamics of 
edge states of the BHZ Hamiltonian when the system is 
suddenly quenched from the TI phase to either the QCP 
or the trivial insulator (TrI) phase. The question here is 
whether there is a surviving edge current following the 
quench. (It is to be noted that when a one-dimensional 
chain of hard core Bosons is quenched from the super- 
fluid to the Mott Insulator state, there is a surviving 
supercurrent in the insulator phase which oscillates in 
time [H]). In our problem, the quench couples the two- 
level subspace of the edge states to a multi-level environ- 
ment of bulk states. We study the decoherence of these 
edge states using the Loschmidt echo (LE) [TS] when the 
2-D Hamiltonian is quenched from the TI phase to the 
QCP (or to the TrI phase). We observe a strong oscilla- 
tion of the low-energy edge states between the two edges 
of the system and the time-averaged persistence of the 



spin Hall current (SHC) when the system is quenched 
to the QCP, which we attribute directly to the linear 
low-energy dispersion at this point (which is also found 
in other models describing 2-D TIs). The experimen- 
tal prospect of real-time tuning of parameters controlling 
QPTs of TIs in optical "50] and photonic lattices [21] as 
well as by exploiting Floquct dynamics [25] , has made the 
study of quenching dynamics of TI Hamiltonians relevant 
and important. It is to be noted that slow quenching re- 
sults in a violation of the Kibble-Zurek scaling of defects 
in systems with edge states [55] . 

The 4x4 BHZ Hamiltonian comprising of two 2x2 
blocks (for opposite electron spins) is given by 



H 



BHZ 



H{k) 
H*{-k) 



(1) 



where H{k) = [C - D{kl + fc2)]l2x2 + A%a'' + kyaV] + 
[m — B{k1 + ky)]a^. Here, A,B,C,D and m are de- 
termined by the thickness of the quantum well and the 
material parameters; the parameter m controls the phase 
of the system and changes sign relative to B when the 
system crosses from the TI phase (where edge states are 
present) to the TrI phase (with no edge states) via a 
DP at m/B = 0. Although the results presented here 
are valid in generic situations, we shall set D = for 
simplicity, which also ensures an electron-hole symmet- 
ric spectrum. We consider a ribbon geometry extending 
from —L/2 to L/2 in the y direction (with the wavefunc- 
tion vanishing at the edges) and apply periodic boundary 
conditions in the x-direction [551 - f57] . 

To obtain the spectrum of Hamiltonian (fl]), we con- 
sider the 2x2 block H{k), fix kx, and use ky — ?► —idy, 
with the trial solution V = (^1,^2)6'*''' [5S]. The con- 
dition that the wavefunctions must vanish at the edges 
of the ribbon quantizes the energies E of the eigenstates 
at a given fc^,, which are given by the solutions of the 
following transcendental equation: 

tanh(A+L/2) tanh(A_L/2) _ 
tanh(A_L/2) ^ tanh(A+L/2) ~ 

Xl + Xl-{B/A)^{Xl-Xl)^ 



A+A_ 
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FIG. 1: Spectrum of the BHZ Hamiltonian in a ribbon geom- 
etry with ribbon width L — 200 nm. The parameters used are 
A = 364.5 meV/nm, B = -686 meV/nrn^ and C = D = 0. 
The values of m used are —10 meV (a), (b) and +10 meV 
(c). There is a small gap of 0{A/L) in the spectrum near 
is = at m = because of the finite width of the ribbon. 
There also exists an exponentially small gap between the edge 
state bands in the TI phase. Note the (almost) equal spacing 
of energy levels for kx = at m, = 0, characterized by the 
solutions of a Dirac particle geometrically confined in a 1-D 
box. 



where A± = -^fe^ +F± yjF"^ - [M^ - E^)/B'^, with 
F = (A"^ - 2MB)/{2B^). The corresponding wavefunc- 
tions are given by 



ipix,y) 



(c+/+(y) + c_/_(y)), 



(3) 



where c± are two-component spinors whose entries are 
determined by the boundary conditions, with f+{y) — 
cosh(2/A+)/cosh(LA-|_/2)— cosh(yA_)/cosh(LA_/2), and 
f-{y) = sinh(yA+)/sinh(LA+/2)— sinh(yA_)/sinh(_LA_ 

/2). 

One can show that in the TI phase [ra/B > 0), there 
are two types of eigenstates of the Hamiltonian, i.e., 
edge states (localized towards the edges and decaying ex- 
ponentially over a length 1/X^(kx,m)) and bulk states 
(spreading across the whole ribbon). In order to have 
true edge states, one needs L ::!> 1/A_. The spectrum, 
which is symmetric in zLk^ and ztE, is displayed in Fig.jT] 
The edge states in the TI phase exist for \kx\ < fcg, fco 



depends on m and L [55] . The solutions of the two 2x2 
blocks are time reversed conjugates of each other, with 
the same set of energies but opposite momentum and 
spin. 

We now perform a sudden quench of the parameter m 
going from m/B > to m/B < and look at the subse- 
quent evolution of an edge state and its spin current. The 
edge states, which originally formed a qubit-like two level 
system with an effective Hamiltonian Hedge ~ Ak^cTz at 
each kx [Ij now get coupled to several bulk modes and 
subsequently decohere. Following a sudden quench, the 
evolution of an edge state is given by 

CO 

\tpedge{kx,t)) = ^ {lpn{kx)\lpedge{kx))e^''^"'^\lpn{kx)) , 

71 — — OO 

(4) 
where \ipedgeikx)) is an edge eigenstate of the Hamil- 
tonian at the initial value of m (= mi) and \ipn{kx)) 
are the eigenstates of the Hamiltonian at the final 
value 7712- The index n runs from — cx) to oo exclud- 
ing n = and denotes the -ve and -l-ve energy bulk 
modes, respectively. Since all the modes are plane waves 
along the x direction, different k^ modes do not cou- 
ple to each other. To study the dynamics of a single 
edge state and quantify its decay, we calculate the LE 

Cit)_= |(V'erfge|e^^(™i)*e-*«(™-)*|V'edge)P, which USiug 



Eq.(|4| can be put in the form 



^edgey'^x^^) 



y^ \{'ipn{kx)\'ipedge{kx))Y 



iE„t 



(5) 

In general, the LE defined above initially drops rapidly 
with time and turns into a rapidly oscillating noisy func- 
tion of small amplitude (Fig. [2la)), indicating that the 
edge state decoheres significantly. However, there is a 
striking difference when one looks at the evolution of a 
low-energy {k^ << k^) edge state following a quench to 
the QCP at m = 0; the LE of the edge state shows a 
pronounced collapse and nearly complete revival for sev- 
eral oscillation cycles (Fig. [2lb)). This is a consequence 
of the nearly equal spacing of the first few energy lev- 
els at low kx near m = (arising due to confinement 
of the linearly dispersing particles (t/ = 1 at the QCP) 
in a ribbon geometry) where the overlap with the edge 
state is the most significant (see Fig. fTl) . We then have 
En ~ sign{n)[Eg/2 + {\n\ — l)AE] for all significant terms 
in Eq.Q, where Eg ^ l/L is the bulk bandgap. The 
summations over n > and ti < 0, then represent Fourier 
series of a peroidic function with period 



27r 

AE 



2% 



A' 



(6) 



making the LE a periodic function with this period. Since 
C{t = Q) — 1, the LE shows a near-complete revival at 
t — riT. Eventually, after several oscillations, the slight 
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FIG. 2: Loschmidt Echo for various edge states and quenches. 
The system parameters are A = 364.5 meV/nm, B = 
-686 meV/nm?, C = D = and L = 400 nm. LE for an edge 
state with: (a) k^ — 0.01 meV/nm and m = — 10 meV after 
quenching to m = +10 meV . There is no significant revival of 
the edge state, (b) k^ ~ 0.001 meV/nm and m = — 10 meV 
after quenching to the QCP at m = 0. There is a pronounced 
collapse to and a nearly full recovery of the LE for several 
cycles, (c) k^ = and m = — 10 meV after quenching to the 
QCP. There is a doubling of the frequency of oscillation as 
compared to the previous case, as this edge state exists on 
both edges. 



non-uniformity in spacing becomes significant and the re- 
vival of the LE weakens. Since AE ^ 1/L, the period 
of this revival scales as L. For small quench amplitudes 
(|m| << A/L), the edge state does not decay signifi- 
cantly. 

Interestingly, we find that the edge state travels from 
one edge to the other and back, existing on opposite 
edges at the points of maxima and minima of the LE 
(Fig. IS^a)). This effect is due to the finite width of the 
ribbon and will not be seen in an infinite system. Since, 
for a significantly large system, the edge states at oppo- 
site edges do not overlap, the LE drops to zero when the 
edge state reaches the opposite edge, and revives again 
when it comes back. For very low values of momentum 
kx — >■ 0, the edge state exists with peaks on both edges 
of a finite ribbon [251. Hence, when a peak on a given 
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FIG. 3: (a) The probability density p — ip^ip of an edge state 
[kx — 0.001 meV/nm) following a sudden quench from m = 
— 10 meV to m — shown at t = (solid) and t = r/2 
(dashed). The edge state travels between the two edges, (b) 
The probabilty current density of the state in the x direction 
(Jo;) at i = (solid) and t — r/2 (dashed). The same state 
carries currents of opposite direction on opposite edges. The 
system parameters are the same as Fig. H. 



edge travels to the opposite edge, the peak on the oppo- 
site edge also travels simultaneously to the given edge, 
resulting in a maximum of the LE at (2?! -|- 1)t/2 instead 
of a minimum, and hence a doubling of the frequency of 
oscillation of the LE. The LE is now minimum at times 
when the state is concenterated near the middle of the 
ribbon (Fig.jl^c)). 
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FIG. 4: The net probabilty current carried by the pair of edge 
states with E — 3.65 meV (near the —L/2 edge over a length 
1/A_ ~ 37 nm) following a quench from m = — 10 meV 
to m — 0. The current collapses and revives like the LE, 
and is always greater than zero, leading to a persistence of 
its time-averaged value. The system parameters are A = 



364.5 meV/nm, B 
200 nm. 



meV/nm , C = D = and L 



The probability current carried near the edge in the x 
direction by the edge states over their decay length 1/A_ 



is proportional to the net SHC carried by the state and 
its time-reversed conjugate in the opposite spin sector. 
It can be calculated using the continuity equation for 
the probability current density J in conjunction with the 
Schrodinger time evolution equation 



{^Hx,y,t)4,ix,y,t))+V-J = 0, 



d_ 

dt 
d 



(7) 



where ip{x,y,t) is the time dependent two- 
component wavefunction of the form ip{x,y,t) = 
{(j)i{y,t),(j)2{y,t))e'^'^^^. Using fc — >• — iV, we obtain 



My,t)^2A{My.tmiy,t) 



(l>*i{y,t)'f>2{y,t)) 



(8) 



with A and B as defined in (fTl). For fc^; — >■ 0, the first 
Dirac-like term dominates the second Schrodinger-like 
term, implying that all low-energy edge states carry vir- 
tually the same current. The profile of Jx is shown at the 
instants of time when the edge state exists on opposite 
edges in Fig.lslb). The sign of the current reverses when 
the state moves to the opposite edge. 

The evolution of the SHC carried by an edge state 
near a given edge (say, —L/2) can be obtained by 
integrating J^ from —L/2 to —L/2 -I- 1/A_; i.e., 

ledgeit) = S-L/2 Jx{y,t)dy. Sincc there are two 

oppositely propagating edge states at a given energy 
which exist on opposite edges of the system, we must add 
the currents carried by both of them. The time evolution 
of such a current following a sudden quench to the QCP 
is shown in Fig. |4] Due to the oscillation of the edge 
states between the two edges, this current also displays a 
pronounced collapse and revival for several cycles. The 
time-averaged value of the current is non-zero and is a 
significant fraction of the original value of the current. 
Thus, there is a persistence of the SHC carried by low- 
energy edge states following a sudden quench to the QCP. 

After several oscillations, the edge state disperses and 
does not regain it's original character again. This dis- 
persion is introduced by the non-linear Bk^ term in the 
Hamiltonian, and must be sufficiently small if sustained 
oscillations are to be observed. Using the condition that 
the spread of the edge state over one time period must be 
significantly smaller than it's initial decay length 1/A_, 
we obtain the condition {BL)/A < 1/A?.. Since we al- 
ready have L ^ 1/A_, this condition implies that these 
oscillations will be seen over an intermediate range of 
ribbon widths roughly determined by the values of the 
parameters A and B. 

It has already been mentioned that recently a realiza- 
tion of a Floquet TI has been experimentally achieved 
in photonic lattices |21j , where the quenching dynamical 



study discussed above can possibly be verified. The dy- 
namics of a Floquet quantum system driven with a time 
period T are determined by its Floquet evolution oper- 
ator U{T), which may be associated with an effective 

Hamihonian H^ff ^ as U{T) = re~^'^^^Jo "^^^'^^ = 
^-[i/h)H^ijT ^ where T denotes the time ordering sym- 
bol. Thus, in a Floquet system, the evolution of an edge 
eigenstate of Hf,ff{mi) following a sudden quench to m2 
is like a normal Schrodinger evolution under Heff{m2) at 
instants of time t — nT. For these oscillations to be vis- 
ible in experimental studies of such a system, one needs 
r<T (see, Eq.Q). 

To summarize, we observe interesting decoherence dy- 
namics of the low-energy edge states of a TI Hamiltonian 
when quenched to the QCP displayed in the temporal 
evolution of their LE. Given the recent prospects of tun- 
ing of control parameters of the related Hamiltonians, we 
believe that this study is of experimental relevance. Fur- 
thermore, one can show that the predicted oscillations 
of edge modes should occur in any 1-D or effectively 1- 
D system with a linearly dispersing [v = 1) QCP (such 
as the 1-D p-wave superconducting chain with Majorana 
boundary modes), due to the fact that the energy levels 
at the QCP in a finite system will be more or less equally 
spaced [29l . 

We acknowledge Apoorva Patel and Diptiman Sen for 
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